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ABSTRACT 

The presence of multiple luminous galaxies in clusters can be explained by the finite time over which a galaxy 
sinks to the center of the cluster and merges with the the central galaxy. The simplest measurable statistic to 
quantify the dynamical age of a system of galaxies is the luminosity (magnitude) gap, which is the difference 
in photometric magnitude between the two most luminous galaxies. We present a simple analytical estimate of 
the luminosity gap distribution in groups and clusters as a function of dark matter halo mass. The luminosity 
gap is used to define "fossil" groups; we expect the fraction of fossil systems to exhibit a strong and model- 
independent trend with mass: ~ 1-3% of massive clusters and ~ 5^-0% of groups should be fossil systems. 
We also show that, on cluster scales, the observed intrinsic scatter in the central galaxy luminosity-halo mass 
relation can be ascribed to dispersion in the merger histories of satellites within the cluster. We compare our 
predictions to the luminosity gap distribution in a sample of 730 clusters in the Sloan Digital Sky Survey C4 
Catalog and find good agreement. This suggests that theoretical excursion set merger probabilities and the 
standard theory of dynamical segregation are valid on cluster scales. 

Subject headings: cosmology: observations — dark matter — galaxies: clusters: general — galaxies: halos 
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1. INTRODUCTION 

Virialized cold dark matter halos grow hierarchically, by 
the merging of smaller virialized halos. Galaxies, which 
populate the halos, also grow hierarchically by the merg- 
ing of pre-existing galaxies. Halo merger rates can be esti- 
mated analytically using the excursion-set theory dBond et all 
which is commonly known as the Extended Press- 
Schechter formalism ( Press & Schechterl ll974t lBoweJfT99 It 
Lacev & ColdfT99l . The merger rates can also be extracted 
from large-scale cosmological simulations. Direct measure- 
ment of merger rates in galaxy groups and clusters, however, 
is challenging. While me rgers can b e identified by unrelaxed 
X-ray morphologies (e.g., Jeltema et al. 2005) or by the pres- 
ence of radio halos produced by nonthermal particles asso- 
ciated with merger shocks (e.g., Enfilin & Rottgering 2002), 
these methods are not yet precise enough to yield accurate 
estimates of merger rates. We here show that a simple statis- 
tic, the luminosity gap, can be applied to large galaxy surveys 
such as the Sloan Digital Sky Survey (SDSS) to test the pre- 
dictions of the excursion-set merger probabilities. 

Dark matter halos merge from the outside in. Mergers ef- 
fectively begin when the components approach within a virial 
radius of each other; they conclude when the separate identi- 
ties of the two halos have been erased, either by the merging 
of the halo centers or by the complete tidal disruption of the 
smaller halo. The duration of the merger defined this way de- 
pends on the rate at which dynamical friction induces orbital 
decay of the tidally truncated subhalos. If a galaxy is located 
at the center of each subhalo, the galaxies appear as separate 
objects until the conclusion of the merger. The timescale of 
orbital decay can exceed the age of the system; this is why 
the number of galaxies in a groups or cluster (the "occupation 
number") exceeds unity. 
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A confirmation of this paradigm can be found in the ex- 
istence of "fossil" groups of galaxies, which contain a sin- 
gle ultraluminous galaxy at the center and no other galaxy 
brighter than L*. The X-ray luminosity of fossil groups is 
comparable to that of rich clusters and indic ates dynamical 
masses - (10 13 -10 14 ) M Q . Uones et alJ 120031) selected fos- 
sil groups with the criterion Amag 12 > 2, where the lumi- 
nosity gap Amag 12 is defined as the difference in photomet- 
ric magnitude between the most and second most luminous 
galaxies in the group; we adopt this definition as well. Fos- 
sil groups have been identified as cluster-sized systems old 
enough that any previous merger with a halo hosting a > L* 
galaxy has completed so that all luminous galaxies have ag- 
glomerated onto the central galaxy. Their incidence rate is 
(8-20)% among observed systems in this mass ra nge, at least 
if the mass is inferred from the X -ray luminosity ( iJones et ail 
2003). Numerical simulations by D'Onghia et al. (2005) pre- 
dict a larger fraction of fossil systems: (33 ± 16)% for mass 
10 14 Mq. These investigations suggest that fossil groups, poor 
clusters, and rich clusters can be distinguished by the time 
elapsed since the last major merger. We show here that the 
incidence rate of fossil groups can be estimated analytically 
from excursion set theory. 

In §|2] we compute the time scale on which dynamical fric- 
tion drives orbital inspiral of a tidally truncated subhalo inside 
the primary halo. In §|3j we estimate the subhalo mass distri- 
bution as a function of the subhalo mass and distance from the 
center of the primary halo. In § |4] we calculate the luminos- 
ity distribution of the most luminous satellite galaxy within 
the cluster and the dependence of the luminosity gap on the 
halo mass. We also derive the fraction of fossil systems in 
groups and clusters as a function of mass. In §|5j we compare 
the predictions of our model to the luminosity gap distribu- 
tionin 730 clusters from the SDSS C4 Catalog ( Mill er et all 
120051) . Finally, in § [6] we show that the observed intrinsic 
scatter in the L c -M relation on cluster scales can be ascribed 
to dispersion in the merger histories. Throughout the paper 
we assume the standard cosmological model consistent with 
the results from the Wilkinson Microwave Anisotropy Probe 
<Snenreleta]J2noa. 
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2. DYNAMICAL EVOLUTION 

Consider a subhalo of mass M s merging into a primary halo 
of mass Mh > M s at redshift z m - We define a "merger" as the 
time when the center of subhalo crosses the virial radius of the 
new composite halo of mass M = Mh +M S . After the halos co- 
alesce, the subhalo spirals toward the center of the composite 
halo. The effective dynamical mass of the subhalo decreases 
after the merger because bound mass is tidally stripped as the 
orbit of the subhalo decays. We denote the bound mass by 
M S (R S ), where R s is the tidal truncation radius. This is related 
to the separation r of the subhalo from the center of the com- 
posite halo via M S (R S )/Rl = M(r)/r 3 , where M(r) is the mass 
of the composite halo contained within radius r. This relation 
yields R s and M S {R S ) in terms of r. 

The subhalo experiences a torque T = |r X F|, where r 
is its position relative to the center o f the composite halo 
and F is force of dynamical friction ( Chandrasekhar [l 9431) 
|F| = 47rG 2 M s (fl s ) 2 p(>-)rn(A)/v 2 . Here, p(r) is the density of 
the composite halo at distance radius r, ln(A) is the Coulomb 
logarithm (which in principle depends on the orbit of the 
subhalo and on the orbital phase space distribution of dark 
matter), and v is the velocity of the subhalo. To simplify 
the calculations we assume circular orbits. Then the veloc- 
ity of the subhalo is the circular velocity in the composite 
halo v = [GM(r) /r\ x l 2 . In numerical simulations of sat ellites 
in hal os ln(A) - 2 (Vela zquez & Whitdll999UFellhauer etafl 
2000). 

The galaxies spiral toward the center of the composite halo 
at the rate 

d JL = J^( d ±Y (!) 
dt M S (R S ) \dr) 

where J = [GM(r)r]^ 2 is the specific angular momentum 
associated with the orbit of the subhalo. We assume 
that the density profile of the halo follows the form of 
iNavarro. Frenk. & White! (1 19971). wit h the concentration pa- 
rameter from the Bullock et al. (2001) fit. Equation Q can be 
integrated to find the distance of the subhalo from the center 
of the composite halo as a function of time. 

3. SUBHALO DISTRIBUTION 

We would like the probability that a subhalo of mass M s is 
located at distance r from the center of the composite halo of 
mass M. The extended Press-Schechter formalism does not 
directly yield such an expression because halos grow through 
an entire hierarchy of mergers. We present a variation of 
one o f the e stablished mode l s for the subhalo mass function 
(e.g.. iNusser & Shethlfl9 99: Fuiitaetal.ll25ol IShethll200l 
lLeell2004t l'O guri & Lee 2004). These models ignore a num- 
ber of issues, such as halo triaxiality and the evolution of sub- 
structure within substructure. Our confidence in their validity 
stems from their success in reproducin g subhalo statistics in 
large-scale numerical simulations (e.g. . IZen tner et aLli2005). 
grav itational lensing observations (e.g., Nataraian & Springel 

2004) , and the cluster luminosity function dCoo rav & Cen 

2005 ) . Note that the known problems with the extended Press- 
Schecht er merger rates are not severe for the mass ratios of 
in terest (tBenson et al.l2005l) . 

Lacev & Cole ( 1993) give an expression for the fraction of 
mass of a halo with mass M at redshift z that lies in progenitors 
with masses between M s and M s + dM s at redshift z m ' 



f(M s ,z m \M,z) : 



(27r) 1 /2[ cr 2( Ms )_ (T 2( M )]3/2 



da 2 (M s ) 



dM s 



x exp 



[Sdz m )-S c (z)] 2 \ 
~2[a 2 (M s )-a 2 (M)]j 



(2) 



where a(M) is the mass variance on scale M and <5 c (z) is the 
critical overdensity for collapse at redshift z- The progenitor 
mass function is obtained by multiplying f(M s ,z m \M,z) by 
t he halo multipli c ity fac tor, dN/dM s = (M/M s )f. 

Lacev & Cole ( 1993) define the formation redshift as the 
time at which the most massive progenitor halo contains half 
of the mass of the final halo. Equation (0 is integrated with 
respect to M s and differentiated with respect to z m to obtain 
the growth rate of the fraction of the halo mass in large ob- 
jects; we interpret this as the probability that a halo at redshift 
z formed at a redshift between z m and z m + dz m , 



dP , N d 

—( Zm ,M,z)= — 

UZ m UZ m J M m 



f(M s , Zm \M,z)dM s (3) 



where M m ; n > lM is a min imum mass cutoff. 

Following Oguri & Lee ( 200J), we multiply the subhalo 
mass function by the formation redshift distribution 

d 2 N dN dP 
~ x (4) 

dM s dz m dM s dz m 

and interpret the result as the probability that the halo acquired 
a subhalo of mass M s at redshift z m - This interpretation is 
only heuristic; doing better requires integrating over merger 
trees. Such a treatment would improve upon the im precise 
definition of "formation time" (Cohn & White 2005) and its 
likely correlation with subhalo properties, but we defer that to 
future work. 

The probability that a subhalo of mass M s lies a distance 
r from the halo center, as a function of the halo mass M and 
redshift of observation z, is given by (we suppress the depen- 
dence on M and z everywhere) 

d 2 N d 2 N dz m 

:7ZT^ M -r)=— — (M s ,zJ^(M s ,z m ) (5) 
dM s dr dMsdzm dr 

for r < R vn -(M,z), where the average density of the halo is 200 
times the mean density of the universe within the virial radius, 
and dzm/dr is the inverse of the derivative of the separation at 
redshift z with respect to the merger redshift z m (from eq. Q, 
after converting cosmic time to redshift). 

The luminosity gap is independent of the subhalo position, 
so we integrate equation with respect to radius to obtain 



dN R 
dM s 



(M s ) = 



Zm(M„0) 



d 2 N 



z m (M, 



S ,R) dM s dz 1 



-(M s ,Z m )dZm, 



(6) 



where R < R vll - is a constant radius and z m (M s ,r) is implicitly 
defined by the relation 



rZmiM.,r) dr 

R vi Az m (M s ,r)]-r= j — dz' m 



(7) 



We solve Q approximately by evaluating the virial radius at z 
rather than z m and by evaluating dr/dz m at R vn -(z)- The latter 
approximation is justified because dr/dz m depends on radius 
only weakly, especially when r ~ /? v ; r . 
The number of subhalos with mass greater than M s is 

n r( m *)= / Sy^VM'. (8) 

Jm s dM' 

Equation (|8j defines a one-to-one relation between M s of a 
particular subhalo and the expected number of subhalos above 
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FIG. 1. — Probability di strib ution for the luminosity gap Amag 12 at z = 
calculated using equation 1141 . From left to right, the curves correspond to 
halos with mass M = (10 1 



10 1 



10' 



2, M m 



i i 



10 15 ' 5 ) M . We assumed In A = 
jM, and R = R VK (M) thi'oughout. While the initial mass of the 
satellite is always smaller than the mass of the primary, the luminosity of 
the satellite galaxy may exc eed t hat of the primary because of the scatter in 
the L c -M relation (equation 1121 . leading to a negative luminosity gap; the 
measured quantity is |Amag 12 |. 



this threshold. It is straightforward to show that the distribu- 
tion of the most massive surviving subhalo mass, M%, is 

if if _..V: ..if, i (9) 



-N R (M 2 ) 



dM 2 dM 2 

Similarly, the distribution of the second most massive surviv- 
ing subhalo mass, M3, is 

dN d 



dM 3 



dM 



■{[l+N R (M 3 )]e- mm) }. 



(10) 



4. THE LUMINOSITY GAP DISTRIBUTION 



We next calculate the distribution of luminosities of the 
first and second most luminous satellites in a galaxy clus- 
ter. Thus we need to relate the initial mass of a subhalo to 
the luminosity L c of its central galax y. At z ~ 0, L c in any 
halo is tightly correlated with its mass (iVale & Ostrikerl2004i 
Coorav & Milosavlievic 2005a b), with a functional form 



L c (M)=Lo- 



(M/M ) a 



'[b+iM/Moyd] 1 / 111 

Our fit to the r-band luminosities in Seli ak et all J2005 ) yields 
L = 5.7x lO 9 L ,Mo = 2x 10 n M Q , a = 4, b = 0.57, c = 3.78, 
and d = 0.23. The relation possesse s lognormal intrinsic scat- 
ter dCoorav & Milosavlievic 2005]J 



dN 
~dL 



(L\M)- 



1 



/27rln(10)I]L 



exp 



log[L/L c (M)] 2 

2S 2 



(12) 



where in the r-band £ « 0.17+JJ of in clusters dCoora\l2 005). 

We assume that the most luminous galaxy is located at the 
center of the composite halo. Then the remaining subhalos 
cannot contain the most luminous galaxy in the cluster; the 
most luminous galaxy in a subhalo is the second most lumi- 
nous member of the cluster. The distribution of galaxy lumi- 
nosities associated with the kth most massive subhalo is then 



dN 
dL k 



dN 



dN 



(L k \M k )^—dM k . 
o dL dM k 



(13) 



The luminosity gap Amag u , is the (observed) magnitude 
difference between the first and the A:th most luminous galax- 
ies in a cluster. The most luminous galaxy is the central 
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FIG. 2. — Probability Pf(M) that a halo of mass M contains a fossil system 
of galaxies. The curves assume In A = 1 and M m j n = jM (solid line), InA = 2 
and M m j n = \M (dashed line), and In A = 1 a ndM m j n = \ M (dot-dashed line). 
The shaded rectangle is the measurement of Jones et al. 1 2003). 



galaxy, and the second one lies in the largest surviving sub- 
halo. In systems of mass M the gap is distributed as 



dN 



d Amag j 



-(AmagjM): 



dN dN , 

—(L k ,M)—(L\M) 
dL k dL 



xS 



Amag lt + -log — 



dL k dL, 



(14) 



where S(x) is the Dirac delta-function. In Figure [2 we plot 
dN /c/Amag I2 for halos of various masses; the median lumi- 
nosity gap is larger in smaller halos. 

The incidence rate Pf{M) of "fossil" systems (Amag 12 > 2) 
can then be calculated by integrating equation (I14> . In Figure 
13 we plot Pf(M) for several combinations of input parameters. 
On scales M ~ 10 I4 M Q , the probability that the sy stem is fos- 
sil is ~ (3-6)%, fairly close to the observed value dJones et all 
2003). It is, however, smaller than the same probability esti- 
mated from the simulations of D'Onghia et al. (2005, see §Q. 
The three curves illustrate the errors we expect from our sim- 
plified treatment; they do not affect our qualitative results but 
may be distinguishable with more extensive observations. 

5. THE LUMINOSITY GAP IN SDSS-C4 CLUSTERS 

We measured the luminosity gap distributi on on 730 clus- 
ters in the SDSS C4 Cluster Catalog dMilleret al J 120051) 
at mean redshift (z) = 0.087. The three brightest clus- 
ter members were identified from the SDSS photometry 
as being within a projected 500/i _1 kpc radius of the cen- 
ter of the cluster. Additionally, these galaxies must have 
m r - nti colors that lie within 2er of the E/S0 cluster ridge- 
line as determined by spectros copically confirmed members 
dVisvanathan & Sandagein"977l) . We utilize extinction cor- 
rected mode l-fit magnitudes and apply z = ^-corrections 
(version 3.2. lBlanton et alJ2003l) . 

The masses of these clusters are estimated from the total 
r-band luminosities via log(/i -I M) « -2.46+ 1.451og(/i~ 2 L r ); 
95% lie in the range (0.5-10) x 10 14 (/z/0.7) _1 M Q . The total 
luminosity is a better mass estimator than the line- of-sight ve- 



locity dispersion or the richness of the cluster ( Miller et al 
2005; see also iLin. Mohr. & Stanford! 120041 iPopesso et al 



2005, and references therein). A histogram of the luminosity 
gap distribution is shown in Figure [3] The mean luminosity 
gap is (Amag 12 ) w 0.75. 

To predict the luminosity gap distribution of the C4 sam- 
ple, we multiply dN /t/Amag 12 of equation dl4> at z = by the 
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Amag 

FIG. 3. — r-band luminosity gap distribution from 730 clusters in the SDSS 
C4 Catalog I Miller et al. 2005). We evaluate the luminosity gap relative to, 
(a) the first and second most luminous galaxies, and (b) the first and the third 
most luminous galaxies. We also show predictions of our model with In A = 1 
and M mm = \M (thick line) and with InA = 1 and M mm = jM (dot-dashed 
line). Fossil systems are located to the right of the dashed line; Pf 2.9% in 
the data and 2.7% in the model. 



mass function of C4 clusters and integrate over mass. The re- 
sulting composite model luminosity gap distribution is shown 
by the thick solid and dot-dashed lines in Figure |3j?. The 
agreement of the model and data is remarkable for the choice 
of the dynamical friction parameter In A = 1. For InA > 2, the 
model distribution develops a local minimum at Amag 12 = 
and its maximum shifts to Amag I2 ~ 0.6 (dot-dashed line in 
Fig. [3}. Such behavior is not apparent in C4 clusters, but it has 
been detect ed in luminous red ga laxies (LRG) in SDSS imag- 
ing data bv lLoh & Str auss ( 2005 ), who measured the luminos- 
ity gap within 1.0/T 1 Mpc of LRGs and found that the peak 
moved to positive values in underdense fields. These curves 
illustrate the sensitivity of the luminosity gap to the detailed 
properties of mergers; clearly the qualitative fit is excellent, 
but more detailed model-fitting may in the future enable tests 
of particular aspects of merger dynamics. 



In Figure[3j> we present the luminosity gap distribution rel- 
ative to the third most luminous galaxy in the cluster (see eq. 
1101 ). The model overpredicts the frequency of small Amag 13 
because it treats the luminosities of the second and third most 
luminous galaxy as independent, whereas the former must ex- 
ceed the latter by definition. 



6. THE ORIGIN OF THE L c -M SCATTER 

The final accretion of a satellite halo at the center of the 
primary halo is accompanied by an instantaneous increase 
of the central galaxy's luminosity. Therefore, the accretion 
history is a source of intrinsic scatter in the L c -M relation. 
We estimate the dispersion E of the central galaxy luminos- 
ity arising through the merger variance alone as the average 
luminosity of the most massive satellite galaxy, (La(M)), as- 
suming no intrinsic scatter in the mass-luminosity relation, 
E(M) - (L 2 (M))/ln(10)L c (M). 1 In our model, 



dN 

(L 2 (M)} = I L c (M 2 )—(M,M 2 )dM 2 



(15) 



We average the dispersion calculated this way over the mass 
function of clusters in our C4 sample; it is £ = 0.23 ±0.01, de- 
pending on the choice of In A. This is not far from our adopted 
value S r = 0.17 and from £ 0.168 measured bv lYang et alJ 
(120031) . Therefore, the dispersion in the L c -M relation in clus- 
ters can be explained by the dispersion in accretion histories. 
If this interpretation is correct, Amag 1/t at fixed cluster mass 
would be correlated with the central galaxy luminosity, im- 
plying that the factors of dN / dL(L\M) in equations dl 3I > and 
dl4> must be replaced by a single multivariate distribution in 
(LjLk); we defer this possibility to a more detailed treatment. 
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1 A more accurate expression, to be presented in a subsequent paper, in- 
volves a sum over the variances of L^M) in the place of (Lo(M)). 
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